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The gravitational clustering hierarchy and dissipative gas processes are both involved in the 
' formation of bulges. Here we present a simple empirical model in which bulge material is 

, assembled via gravitational accretion of the visible companion galaxies. Assuming that merging 

' leads to a starburst, we show that the resulting winds can be strong enough that they self- 

regulate the accretion. A quasi-equilibrium accretion process naturally leads to the Kormendy 
t I . relation between bulge density and size. Whether or not the winds are sufficiently strong and 

' long lived to create the quasi-equilibrium must be tested with observations. To illustrate the 

model we use it to predict representative parameter dependent star formation histories. We 
find that bulge building activity peaks around redshift two, with tails to both higher and lower 
' redshifts. 

(N 

^ ■ 1. Introduction 



Bulges are stellar dynamical pressure supported systems that generally have much 
higher surface brightnesses than galactic disks. They therefore have undergone more 

■ collapse than galactic disks, evidently with the angular momentum barrier removed. 
' Galaxy merging is an inevitable process that redistributes any pre-merger stars into a 

physically dense, but phase density lowered, pressure supported distribution. Stellar 
dynamical mergers produce objects with flattenings largely unrelated to their rotation. 
• In the presence of gas, merging is empirically associated with an often dramatic rise in 

star formation. These new stars that are formed in place almost certainly reflect the 
Q ■ chemical history and the dynamical state of the growing bulge. This paper calculates 

| H I some of the properties of bulges expected on the basis of merging with star formation of 

c/3 ■ largely gaseous pre-galactic fragments. 

The rate of major mergers can be calculated directly from the observed numbers of 
close pairs of galaxies. Remarkably, this is now an observational quantity for which there 
are measures from low redshift to the "U-dropout" population centred around redshift 

■ three. There are some substantial uncertainties in the various observational quantities 
I which go into the merger calculation. The details of this calculation will become much 

more precise over the next few years as the evolution of the two-point correlation function 
becomes better determined. 

Mergers are widely observed to induce an intense nuclear starburst. Theoretically 



this is at least partially understood ( Barnes fc Ilernquist 1992 ) on the basis that the 



strong dynamical interactions during a merger leads to a loss of angular momentum in 
a cool gas, helping to funnel it to build up a dense central gas reservoir from which 
stars form at astonishing rates in a starburst. Starbursts in turn develop winds which 
we suggest can lead to the accretion being a self-regulating, although this is dependent 
on the ram pressure and duration of the wind. Moreover, self-regulating accretion can 
lead to quasi-equilibrium star formation in the bulge which can lead to some of the 
observed regularities of bulge properties with mass or size. At this stage the details 
of this picture are speculative, but are open to observational refinement, which helps 
motivate the calculations presented here. 

Other papers presented at this meeting describe in detail the properties of bulges. 
Here we take the properties of "classical bulges" to be roughly as follows (e.g. 



Wyse Gilmore fc Franx 1997) 
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• Bulges follow the "Kormendy relations" , that is, the characteristic surface 
brightness correlate strongly with the scale radius (Kormendy 1977, De Jong 1996, 



Pahre, Djorgovski & De Carvalho 1995). 



• The flattening of the figure of bulges is approximately consistent with their rotation 
( pavies et al. 1983D . 

• Bulge stellar populations are predominantly old, although there are well-documented 
cases of relatively young bulges. 

• Bulges have a mass-metallicity correlation. 

A useful model for bulges should be able to provide a physical origin for these properties. 

The paper has three main sections. In §2 we discuss the empirically determined rate at 
which mergers occur as a function of redshift. Then, in §3, we discuss star burst winds, 
and the effects those winds will have on accreting gas. Section 4 pulls these two together 
in some specific model calculations. 



2. Merger Rate Measures 

A host galaxy has a number of near neighbors within radius r and pairwise velocity 
\v\ far above the mean density no (in proper co-ordinates), 

N{<r,<v)^ATTno{l + zf / £_{r\z)f{v\z)r^v^ dr dv, (2.1) 

Jo Jo 

where f(r|z) and f{v\z) are the redshift dependent two point correlation and velocity 
distribution functions, respectively. We have made the important assumption that the 
distribution of pairwise velocities is constant over the separations of interest. This is 
not true in general, but is sufficient for our application to the relatively small scales, 
r < 20 kpc, that are of interest for merging. To calculate merger rates we need estimates 
of the correlation function on small scales, the pairwise velocity dispersions, and the 
mean time for a merger to occur within this volume. 

2.1. Close Pairs and the Correlation Function 

The galaxy-galaxy correlation function is accurately modeled as a power law, ^(r) = 
(jo/fV- The reliability of this power law on scales less tha n about 100/t~^ kpc 



relevant to galaxy mer ging is discussed in detail for the SSRS2 ( Patton et al. 1999a ) 



and CN0C2 samples ( |Patton et al. 1999bD . These papers support three important 



conclusions. First, the power law extrapolation of the correlation function to 2Qh^^ kpc 
is consistent with the density of pairs measured inside this radius. Second, the R 
band luminosity function of galaxies in 20/i~^ kpc pairs is consistent with being drawn 
from the field luminosity function. This property allows the fully general correlation 
function, ^(r, W|| , uj^, Li, L2), to be factored into a luminosity function, and a spatial 
and kinematic correlation function, (j){Li)(j){L2)(,{r)f{v) where we drop the distinction 
between velocities along the line of separation, and perpendicular velocities, t;j_ 



(Peebles 1980). The luminosity factorization glosses over the various lines of evidence 



(Carlberg et al. 1998, Loveday et al. 1995) that the theoretically expected weak increase 
of correlations with galaxy mass does exist in the correlations. However, this relatively 
small effect cannot yet be detected in the current samples of close pairs which have not 
yet broken through the barrier of 100 pairs. The third important result is that there 
is morphological evidence that r < 20h~^ kpc pairs are indeed interacting at a level 
that indicates that these are high probability mergers-to-come. The 20h~^ kpc scale is 
also chosen such that the galaxies are not so strongly interacting that their luminosities, 
morphologies and colors bear little resemblance to their unperturbed values. 
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Figure 1. Measured co-moving correlation lengths as a function of redshift. The points are 
the LCRS at low redshift, CN0C2 at intermediate redshift, and the Giavalisco et al. (1998) 
close pairs and the "box counts" of Adelberger et al. (1998) at redshift 3. The li nes are for 



the correlation lengths measured in simulations ( Colin, Carlberg fc Couchman 1997 ) for galaxy 
halos (solid lines) and the mass field (dashed), with = 0.2, and Qa ~ (heavy) or Qa = 0.8 
(light). 



The volume integral of the power law correlation function, ^(r|z) = (ro(z)/r)''', is, 

An ( ro(2:) 



47r 



^(r|z)r dr — 



3-7 



(2.2) 



The redshift dependence of the average density inside a. r — 2Qh ^ kpc neighborhood 
around a galaxy is estimated using the preliminary CN0C2 correlation 7 — 1.8, 

ro(z) = 5.15(1 + z)[^"'^+')/''l/i"^ Mpc (co-moving), (2.3) 

where e = —0.6 ± 0.4 ( [Carlberg et al. 1998 ). Using this in Equation ^.2| we find that 
the integrated density inside 20/i~^kpc is 1.56[(1 -I- z)/1.3]~'^no Mpc (proper units). 
In Figure |^ we show the co- moving correlation length as a function of redshift. We 
also plot the correlations of the galaxy mass halos and the particles in simulations 
( polin, Carlberg fc Couchman 1997 ). It should be noted that there is a relatively good 
understanding of why the halo correlation function shows relatively little evolution 
( parlberg 1991 , Baugh, et al. 1998). Observationally it is currently acceptable to take 
Tq to be fixed over this redshift interval, or, e = —1.2 for a 7 = 1.8 power law. 

2.2. Pairwise Velocities 

The CN0C2 velocity distribution function, f{v), at small scales is acceptably 
modeled as a Gaussian, although an exponential also provides a similar quality fit 
( parlberg et al. 1998 ). We will take the velocity to be isotropic, 0-3 = Vicr^, where 
az is the velocity dispersion measured along the line of sight. Redshift surveys can 
be used to measure (Tz. At separations of about 1 Mpc the the pairwise velocity 



dispersion in the redshift direction is about 300 km s ^ (e.g. Davis & Peebles 1983, 
Marzke, Geller, Da Costa & Huchra 1995, Carlberg et al. 1998), constant with redshift 



over the z < 0.5 range. If the critical velocity to merge is taken as about 220-\/2 



4 



R. G. Carlberg: Bulge Building with Mergers and Winds 



kms^^ then the fraction of close pairs with sufficiently low velocities to merge is 
merely 5.1%. This does not accord well with the impression that most close pairs 
have such large tidal tails that they are almost certainly doomed-to-merge pairs 



(Toomre & Toomre 1972, Tooomre 1977). In a non-merging, equilibrium distribution 
the pairwise velocity declines as oc r^"'', where 7 is the slope of the power law 
correlation. The one dimensional pairwise velocity dispersion at 20 kpc is therefore about 
190 km s-^ 

The dynamical details of pair mergers in a cosmological setting that includes 
the tidal fields of surround ing structure have not been studied in detail (but see 



Carlberg fc Couchman 1989 ), such that one of the best estimates of the critical velocity 
to merge remains Aarseth and Fall's (1980) value of Vmg = 1.2-\/2up, where Vp is the 
"parabolic" velocity, Vp = fpVc, at the orbital pericenter, where fp is at least \/2 (for a 
point mass). This "parabolic" velocity is that for those pairs assured to merge, in the 



absence of tidal fields ( fToomre fc Toomre 197^ [Aarseth fc Fall 1980[ ). For /p = \/2 the 



fraction of the the close pairs with velocities low enough to merge is 27% but this rises 
quickly with increasing /p, becoming 41% and 54% at fp = V3 and 2, respectively. We 
will use an fmg — 0.5 but recognize that this number is both empirically and theoretically 
uncertain. 

2.3. Merger Times of Close Pairs 

As a reference timescale for merging we start from our reference radius of 20 kpc where 
the time for a circular orbit is 0.62 Gyr at a speed of 220kms~^. In detail the rate 
of inflow through the 20 kpc depends on mass and the orbital details, so we use a 



merger timescale at our reference radius of 20 kpc of 0.3 Gyr (Barnes fc Hernquist 1992 



Dubinski, Mihos fc Hernquist 1999). We use our estimate that fmg = 0.5 of all 20 kpc 
pairs have pairwise velocities less than the critical velocity for merging. If the pairwise 
velocities and substantially higher than the critical velocity for merging then the merging 
fraction drops nearly as v^, which is such a dramatic change that it should be testable 
via the morphology of galaxies in close pairs. 

2.4. Estimated Merger Rates 

Combining our estimates of clustering, pairwise velocities and the available understanding 
of the dynamics of merging, we find that the specific merger rate is 2.4[(1 + 
z)/1.3]-^noMpc^Gyr"\ If we take uq as being the CN0C2 galaxy number density 
to 0.21L* adjusted to O.IL,, uq ~ 1.1 x 10"'^ Mpc~^ (co-moving) the merging event rate 
for galaxies above our minimum mass is 1.3 x 10^^(1 + z)^'^ Gyr^^. The rate of accretion 
of pre-merger stars onto galaxies as a result of merging is relatively slow. The time scale is 
L*/7?,L ~ 60 Gyr, or about 5 Hubble times at z = 0.3. This result is based on the directly 
visible galaxies, L > 0.21L,, which contain about 80% of the total stellar luminosity 
galaxies. This relatively low rate of accretion of visible galaxies relative to the Hubble 



times continues on to the z ~ 3 regime (Giavalisco et al. 199S, Adelberger et al. 1998), 



where the number densities and the co-moving correlations are similar to those observed 
for present day galaxies. That is, for go = 0.1 their data indicate a density of 7?, < 25.5 
mag Lyman break galaxies of uq ~ 2.2 x lO^'^Mpc"'^ and a co-moving correlation length 
of 5h~^ Mpc. The self-event rate of this population is 2.0 x~^ Gyr~^, only 3.8 times that 
at z ~ 0.3. Since the cosmic time at z ~ 3 is about 20% of that at z = 0.3 the relative 
impact on the hosts of these self-mergers is small adding perhaps 5-10% more mass over 
the entire z = to 3 range. For the lower luminosity galaxies (about 2 magnitudes 
fainter than those with spectroscopic redshifts) inferred to be in this redshift range 
from the HDF the volume density is about a factor of 20 higher, but the implied cross- 
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correlation length is about a factor is about a factor of \/3 smaller (Bteidel et al. 1999), 
where we assume that the cross-correlation depends on the product of the relative biases. 
This implies that the high-low luminosity merger rate is about 0.5 Gyr"^ which is large 
enough to build a galaxy over a z ~ 1 — 4 interval. In the intermediate redshift range 
the M > 0.2il/* galaxies cannot self- merge to significantly alter the mass function. At 
higher redshift the lower luminosity Lyman break galaxies rise very steeply in number, 
a ~ —1.8 ( ptcidcl et al. 1999| ). These large numbers completely change the situation, 
allowing their mergers to substantially alter galaxy masses. 



3. Starburst Winds 

There is a remarkably simple physical description of what happens when star formation 
is rapid in a relatively small volume. The inevitable outcome is a very strong wind. 
Chevalier & Clegg (1985, hereafter CC) simplify the situation to the equilibrium solution 
of mass injection at a rate M, with accompanying energy injection, E, in a sphere of 
radius R. Since it turns out that the wind velocities exceed 1000 km s~^, gravity can be 
ignored in a first approximation. CC provide a full solution at all radii. Here we are 
mainly interested in the asymptotic solution at large radii, where we cast the CC solution 
in terms of the terminal wind velocity, u and the mass injection rate, M^, 

This wind produces a ram pressure of, 

P^=p^u (3.5) 

where R is the size of the region into which the mass and energy are injected. 

As representative numbers we will take R = 10^^ cm, about 1/3 kpc, and = 
lO^^gs"^, about 15MQyr~"'^, approximately the mass injection rate expected during a 
burst of star formation of 1500 Moyr"'^. We follow CC and use u = 2000 kms~^. In the 
central region, 

~ 0.296^, (3.6) 

or 1.5 X lO^^^'g, or a number density of about Icm^'^. 

The cooling time, tcooi = SkT/{nA{T)), at T — lO^K where the cooling rate is about 
A ~ 3 X 10~^^cm'^s~^ is about 4 x 10^ yr. The flow time across the region is only 
1.6 X 10^ yr, so the hot wind does not have time to cool. Denser regions in pressure 
equilibrium will cool more quickly so that the ISM is unstable and bound to consist of 
the hot wind phase and one or more cool phases. Many of the aspects of this situation 
are discussed in Ikeuchi & Norman (1991). 

3.1. Ram Pressure Stripping 

A major objective of this paper is to estimate the ram pressure stripping of the hot 
wind on infalling objects. We calculate the effects of ram pressure stripping, but note 
that transport processes, such as turbulence and heating, could help to increase the rate 



of gas removal (Nulsen 1982). The calculation proceeds in a series of steps. First, we 
derive for our specific case the fairly standard results that the wind will have a very 
high momentum flux. The wind has a sufficiently low density that it will move out 
before cooling. The infalling objects are taken to be angular momentum conserving, 
but maximally dissipated disks in dark halos, approximated as truncated Mestel disks. 
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The fractional stripping can be easily estimated for these objects. The strength of the 
starburst wind is calculated under the assumption that star formation is occurring on 
timescale comparable to the crossing time of the bulge. Bringing these elements together 
gives an expression for the fractional mass of an infalling object which succeeds in joining 



the bulge, Equation 3.16 



The ridge line of Kormendy's relation is = 3.021ogro + 19.74 B mag arcsec 



-2 



(Kormendy 1977, De Jong 1996), where /^s is the B band surface brightness and ro 
is the bulge scale radius in kpc. For a constant mass-to-light ratio this translates to 
p cx r^'^''^ . This density radius relation would be slightly weakened if we allowed for a 
decrease in mass-to-light for lower luminosity systems. The implied densities are high 
enough that bulges are self-gravitating. If the rate of infall is, on the average, regulated 
by the starburst winds, then the fact that the infall rate is very insensitive to the mass of 
the host implies that the characteristic radius and density of the bulge will scale roughly 



as p cx R^^ , as in Equation 3.14 . This relation is roughly the ridge line of the Kormendy 
relation. This is physically easy to understand. The total bulge star formation rises as 
the bulge gas density, ph, to the 3/2 power. For a given accretion rate, a rise in bulge 
density will increase the SFR, and hence the outgoing wind, which temporarily reduces 
accretion, allowing the gas density to be reduced. 

The ram pressure rises as r~^ with decreasing radius. At i?, the outer radius of the star 
forming volume, the surface density below which stripping occurs is rises to its maximum, 

4^ SI' 

where Vc is the circular velocity of the incoming gas in its dark halo. For the starburst 
numbers above, E = 0.16 gcm~^ is the maximum surface density that can be blown away 
via ram pressure alone. For comparison, the central surface density of a disk galaxy is 
about lgcm~^. We conclude that the effects of this wind would be significant even on 
disk like that of the Milky Way if it were completely gaseous. 

If we collapse the gas inside an isothermal halo with a velocity dispersion as and an 
angular momentum parameter A ~ 0.05 to centrifugal equilibrium its surface density 
increases by about a factor of 10^ over that of the projected isothermal halo, 

10'f^6^, (3.8) 

where d is the radial co-ordinate in the disk. For our typical case, Q.b — 0.01, 
H = 100 km s""'^ Mpc^^ we find that this Mestel disk has a surface density profile. 

The total mass in a Mestel disk of the above form diverges. If we cut the disk off at the 
radius of the last orbit that can have come from the outermost virialized part of the halo 
we can estimate a total mass. The halo is virialized inside approximately r200j the inside 
of which the mean density is 200/5c(-2), where Pc{z) = 3H{z)'^ / (SttG). For an isothermal 
sphere, M{r) — 800/37rpc{z)r2QQr. The isotropic velocity dispersion that maintains this 
sphere in equilibrium is, = 400/37rGpc(2)f'2oo- With the definition of Pc{z) we find 
that for the isothermal sphere. 

This is a physical (proper) distance. In this potential the circular velocity is = y/2a. 
The total gas mass in the disk (assuming that there are few stars) is 200pcf^6|7T2oo, 
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or, 



which is needed to find as which is required in the evaluation of d/r for stripping below, 



Equation 3.16 



3.2. Starhurst Rates 

Both empirical evidence and theoretical considerations suggest that the timescale 
for star formation should be proportional to the local dynamical timescale 
(Lchncrt & Hcckman 1996, Kennicutt 199^ ). The available data point to a relation of 



tsfr — R/vb, where vt is the local circular velocity, which in the case of a stellar bulge may 
be due to self-gravity, not the dark matter background. Therefore Msfr = Mvt/R, and 
we will take the wind as — ewMsFR- The mass can be expressed as M = An/SpR^, 
where p refers to the gas density. 

The wind that results from this starburst has a ram pressure, 

47r 

Pw = ewVbup—^. (3.12) 
3 

3.3. Self-regulating Starbursts 
The starburst wind will blow away infalling surface densities smaller than, 

47r i?2 

T,— < e^Vbup— — . (3.13) 
r 6 r'^ 

In equilibrium this leads to a balance of the bulge which creates the wind and the and 
infall, 

pR^^^^_Ytm_ (3.14) 

AtTEwU Vb 

If bulges are to be self-gravitating, then Vb > Vc, and generally they are found to have 
circular velocities comparable to those of the disk. The important thing to note is that the 
quantity pR^ is completely determined by the starburst, whereas the RHS is completely 
determined by the physics of infall. 

For Mestel disks derived from satellite halos of Cs, and = 2ttGpR^ are stripped at 
radii beyond, 

- > 47^-10^^^., (3.15) 

r 4e^ u 

For typical numbers we take as = ah/2 and e^, ~ 10^^. We find stripping at d/r > 1.3, 
which for r = R, and i? ~ 1/3 kpc is 99.5% of the disk mass. 

The equation for stripping of Mestel disks can be used to predict the rate of successful 
accretion in the presence of a starburst. We simply multiply with the ratio of d/r 
for stripping. The accreted mass is simply 

d R 

Maccrctc = -77; Mm- (3.16) 



4. Realizations of Bulge Formation Histories 

Combining merging rates and stripped fractions to build bulges with redshift is now 
straightforward. To illustrate the model we do a simple Monte Carlo simulation. We 
start the simulation at when the universe is about 0.5 Myr old, prior to significant galaxy 




O I ' 1 ' 1 ' 1 ' 1 ' 1 

10 8 6 4 2 
redshift 

Figure 2. Mass redshift history for 100 realizations of bulge building. The calculation assumes 
that suitable pre-galactic objects, disks in dark halos, are available at the beginning of the 

calculation. 




redshift 



Figure 3. Fraction of infalling gas mass successfully accreted versus redshift for the 100 

realizations of bulge building. 



creation or merging. We adopt an Hq = 65, Q,m = 0.2, A = cosmology, although the 
results are not very sensitive to the precise choice of cosmology. The mergers occur at a 
rate independent of mass at at a rate 7?.„(1 + Time is divided into 0.5 Myr intervals 
and the probability of an accretion event is the rate per unit time multiplied with the 
time interval. An accretion event occurs if a [0,1] random number generator produces a 
number less than this probability. The accreted objects begin with a mass drawn from 
4>{M) oc M"e~^ , where the normalizing constants are unity. This therefore assumes that 
the Mestel disks in their dark halos are largely present when the calculation is turned on. 
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Table 1. Merger-Wind Simulation Parameters 



Parameter 


Symbol 


Default Value 


Minimum satellite mass 




0.1 


Mass function slope 




-1.8 


Current merger rate 




0.01 


Merger-redshift index 


m 


1.2 


Gas fraction 


^gas 


0.01 


Bulge size 


Rb 


10^^ cm 


Wind duration 


T 


10* yr 


Mass loss efficiency 




0.01 
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o 
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o 
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00 



o 

1- o 
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o 
o 



o 
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50% moss redshift 
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redshift 

Figure 4. Redshift of 50% mass assembly for 10000 realizations of the standard model. 



Clearly this is not accurate at larg e redshift, but is ar guably a useful assumption over the 
redshift zero to about four range ( ^teidel et al. 19991 ). There is no presumption that the 
characteristic mass is as large as that of a full galaxy today, however the characteristic 
mass does need to be comparable to a bulge mass, since we find that the final masses 
are distributed around unity, the characteristic mass of the infalling objects. The wind 
is assumed to blow uniformly for a duration of from the previous merger. If a new 
accretion event is generated during this interval, then the mass of the incoming satellite 
is reduced according to the stripping equation. The parameters and their default values 
are outlined in the following table. 

In detail there are many (mostly minor) complications that are swept under the rug 
here. The model is quite naive in that we assume that there is a ready supply of gas 
containing companion galaxies with roughly a galaxy mass distribution. Observation 
seems to support this as being true from redshift zero to about four, which covers most 
of the activity here. Mergers wouldn't do much at all if a > —1-2, as is observed at low 
redshift. We have chosen a steeper a to both mimic the increase in gas content with 
decreasing mass, and, to take into account that a does become steep in the redshift 3-4 
range (Steidel et al. 1999). The gas content of galaxies should likely vary with redshift, 
whereas we have simply taken them to be all gas. Likely the pairwise velocity dispersion 
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redshift 



Figure 5. Redshift of 20% mass assembly for 10000 realizations of the standard model. 
Comparison with Figure ^ emphasizes the very large spread in redshifts of formation. 



decreases somewhat with increasing redshift, which will diminish the fraction of low 
redshift pairs that merge. Overall the effect of all of these things would be to decrease 
the importance of merging at below redshift one. The natural tendency of the model is 
to have little activity at low redshift anyway, so the basic character of the results should 
not depend on these simplifications. In any case, the main purpose of this "toy" star 
formation history is to examine the basic viability of the model, not to fine tune the 
parameters. 

Overall, we find several interesting results. First, the mass buildup predicted by this 
simple model seems to be very roughly in accord with the requirements of "classical" 
bulges. Figure ^ shows that About 10% are more than 50% formed at redshift 5. But, 
about 3% are only half formed at redshift 0.5. Although the median time of half assembly 
is a reassuring redshift of about two, there is a tremendous spread of formation times. 
The redshifts of 20% assembly are shown in Figure |[ 

The standard wind has a lifetime of 10^ years, which in many cases limits infall to about 
2/3 of what it would normally be. Figure ^. The limitation of infall would help drive 
bulges towards the Kormendy relations. Winds are even more effective if the smallest 
mass to be merged into a bulge is reduced to 3% of M,, rather than our standard 0.1M». 
This is shown in Figure ^. The basic merger rate is the same, so the redshifts of assembly 
are not greatly altered. However because the bulge is being built of more, smaller, units, 
the buildup has less dispersion in time. 



5. Conclusions 

The bulge formation history is predicted here using the observed density of nearby (gas- 
rich) galaxies with masses comparable to the eventual bulges. Further, we have argued 
that starburst winds will have a significant effect on the accreting gas. A straightforward 
assessment of our results is that the merger history appears to be roughly in accord with 
what is known about star formation histories and bulge ages. An additional step, not 
taken here, is to use these assembly histories to predict the color distribution of bulges as 
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redshift redshift 



Figure 6. The assembly of a bulge built from pre-galactic objects with a minimum mass of 
O.OSAf, , rather than the standard 0.1i\f, and wind lifetime of 10^ yr. The left panel is the analog 
of Figure ^ and the right panel is the analog of Figure |[ The assembly here is smoother and 
more affected by the starburst winds. Although there is less dispersion in formation times, the 
median 20% and 50% time are not significantly different from the standard model. 



a function of redshift. These then become a simple but powerful test of the model. The 
attraction of the merger model is that it is based on observations, which now provide 
a very basic measure over the redshift zero to four range, although the masses of the 
incoming objects are not well quantified at the highest redshifts. Beyond redshift four, 
this model is likely becoming less reliable, since the assumption that the pregalactic 
fragments are in place has no basis in observation and even in a low density Universe, 
the Press-Schechter approach would indicate that the halos in which these reside are 
becoming less numerous. Furthermore the gas content probably evolves with redshift. 
An attraction of this approach is that the need for additional physical parameters can 
be driven by observation. 

The relevance for bulges of the strong starburst winds that we have advocated is far less 
clear at this stage. The attraction of the idea is that it promotes the development of the 
Kormendy relations in a self-gravitating system. The tests of this model will require fairly 
detailed observations of bulges at high redshifts, particularly concentrating on the heavily 
reddened objects in which starbursts generally occur. The presence of starburst winds 
and their affects can potentially be detected in optical emission lines. X-ray observations 
of hot gas (e.g. Heckman et al. 1999[ ) and ultimately in resolved observations of the 



molecular gas (e.g. Prayer et al. 1999| ) which is known to exist. 



This research was supported by NSERC of Canada. I thank the Carnegie 
Observatories, Pasadena, for their hospitality during the time when this work was 
initiated. 



REFERENCES 
Aarseth, S. J. & Fall, S. M. 1980, ApJ, 236, 43 

Adelberger, K. L., Steidel, C. C, Giavalisco, M., Dickinson, M., Pettini, M., & Kellogg, M. 
1998, ApJ, 505, 18 

Barnes, J. E. & Hernquist, L. 1992, ARAA, 30, 705 



12 R. G. Carlberg: Bulge Building with Mergers and Winds 

Baug h. C. M. Benson A. J., Cole, S., Frank, C. S., & Lacey, C. G. 1998, MNRAS, submitted 
( |astro-ph/9811222j ) 

Carlberg, R. G. & Couchman, H. M. R 1989, ApJ, 340, 47 

Carlberg, R. G. 1991, ApJ, 367, 385 

Carlberg, R. G., Yee, H. K. C, Morris, S. L., Lin, H., Sawicki, M., Wirth, G., Patton, D., 
Shepherd, C. W., Ellingson, E., Schade, D., Pritchet, G. J., & Hartwick, F. D. A. 1998, 
Phil. Trans. Roy. Soc, in press. 

Chevalier, R. A. & Clegg, A. W. 1985, Nature, 317, 44 

Colin, P., Carlberg, R. G., & Couchman, H. M. P. 1997, ApJ, 390, 1 

Davies, R. L., Efstathiou, G., Fall, S. M., lUingworth, G. & Schechter, P. L. 1983, ApJ, 266, 41 
Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 

De Jong, R. S. 1996, A&Ap, 313, 45 

Dubinski, J., Mihos, J. C, & Hernquist, L. 1999, ApJ submitted ( |astro-ph/9902217| ) 
Prayer, D. T., et al. 1999, ApJL, 514, L13 

Giavalisco, M., Steidel, C. C, Adelberger, K. L., Dickinson, M. E., Pettini, M., & Kellogg, M. 

1998, ApJ, 503, 543 

Heckman, T. M., Armus, L., Weaver, K. A., & Wang, J. 1999, ApJ, in press ( [a.stro-ph/9812317| ) 
Ikeuchi, S. & Norman, C. A. 1991, ApJ, 375, 479 
Kennicutt, R. C. , Jr. 1998, ApJ, 498, 541 
Kormendy, J. 1977, ApJ, 218, 333 

Lehnert, M. D. & Heckman, T. M. 1996, ApJ, 472, 546 

Loveday, J., Maddox, S. J., Efstathiou, G. & Peterson, B. A. 1995, ApJ, 442, 457 
Marzke, R. O., Geller, M. J., Da Costa, L. N. & Huchra, J. P. 1995, AJ, 110, 477 
Nulsen, P. E. J. 1982, MNRAS, 198, 1007 

Pahre, M. A., Djorgovski, S. G. & De Carvalho, R. R. 1995, ApJL, 453, L17 

Patton, D. R., Marzke, R. O., Carlberg, R. G., Pritchet, C. J., & da Costa, L. N. 1999, ApJ, to 
be submitted. 

Patton, D. R., Carlberg, R. G., Marzke, R. O., Yee, H. K. G., Lin, H., Morris, S. L., Sawicki, 
M., Wirth, G. D., Shepherd, C. W., EUingson, E., Schade, D., & Pritchet, C. J., 1999, 
ApJ, to be submitted. 

Peebles, P. J. E. 1980, Large Scale Structure of the Universe (Princeton University Press: 
Princeton) 

Steidel, C. C. Adelberger. K . L., Giavalisco, M., Dickinson, M, & Pettini, M. 1999, ApJ, in 
press ( [a,stro-ph/981139E| ) 

Toomre, A. & Toomre, J. 1972, ApJ, 178, 623 

Toomre, A. 1977, in Evolution of Galaxies and Stellar Populations, ed. B. M. Tinsley & R. B. 
Larson, (Yale Observatory: New Haven) p. 401 

Wyse, R. F. G., Gilmore, G. & Franx, M. 1997, ARAA, 35, 637 



